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1 . Introduction 

The approach to the development of MODEL III has been to 
divide the problem into three steps of increasing complexity. In 
Chapter IV we successfully developed a variational algorithm for 
the classical temperature retrieval problem that includes just the 
radiative transfer equation as a constraint. The radiances for each 
of the four TOVS MSU microwave channels were dependent variables. 

Chapter V summarizes the second step which combines the four 
radiative transfer equations of the first step with the equations 
for a geostrophic and hydrostatic atmosphere. This step is intended 
to bring radiance into a three-dimensional balance with wind, 
height, and temperature. The use of the geostrophic approximation 
in place of the full set of primitive equations allows for an 
easier evaluation of how the inclusion of the radiative transfer 
equation increases the complexity of the variational equations. 

It should be noted that the variational method is a powerful 
mathematical tool and a powerful method for diagnosing the physical 
role of the observations in the adjustment. We developed seven 
different variational formulations for the geostrophic, hydrostatic 
and radiative transfer equations. The first derivation was too 
complex to yield solutions that were physically meaningful. For 
the remaining six derivations, the variational method gave the same 
physical interpretation - the observed brightness temperatures 
could provide no meaningful input into a geostrophic, hydrostatic 
balance - at least through the problem-solving methodology employed 
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in these studies. It would be axiomatic therefore, that the 
brightness temperatures could provide no meaningful input into a 
variational assimilation with the primitive equations. 

During the writing of this chapter, the equations were 
reviewed and a conceptual error regarding one of the Lagrange 
multipliers was discovered. 

In the following section, the variational methodology is 
presented and the Euler-Lagrange equations rederived for the 
geostrophic, hydrostatic and radiative transfer equations. Then 
the equations are reduced in number through elimination of 
variables to produce a single equation for the geopotential height. 
It is shown that the single equation is too difficult to solve but 
that a three equation set can be solved iteratively. It is also 
shown that space-based thermodynamic data can be assimilated into 
the meteorological data mainstream and that none of the 
difficulties associated with traditional temperature retrievals 
will be encountered. 

2. A Variational Assimilation Theory for the Geostrophic, 

Hydrostatic and Radiative Transfer Equations 

The variational formalism will be derived for the four 
radiative transfer equations in integral form. Let the dynamical 
constraints be, 


^-(JW^+p-o 


( 2 ) 
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m 3 -v-<b x -F s -0 (3) 

/n 4 =u+<J> y +F 6 =0 (4) 

For additional simplification, set the terrain correction term 6=0. 
The forcing functions F 5 and F 6 (see Chapter II) are simplified 
through setting R 0 = 0. 

The integrand of the functional to be minimized is, 

I-k-l ( u-u °) 2 +ti 1 (v-v°) 2 +tc 2 (T-T°) 2 +ti 3 (4>-<t>°) 2 
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where the n i are the relative weights accorded to the observations. 

Performing the variations for the eight dependent variables, 
u, v, §, T, and Bj (j=l,4), yields the following Euler-Lagrange 
equations, 


6 u: n x (u-u °) +X 4 -0 


( 6 ) 
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dv: n 1 (v-v°) + A 3 -0 (7) 

6<(>: n 3 (4>-<J>°) -X 2a +X 3x -X 4y -0 (8) 

6r: % 2 (T-T o )+Yk 2 -^X 1J J^w j da-0 (9) 

85: g [n 4J (flj-fl/)+X w -0] (10) 


These eight equations plus the seven original constraints 
constitute a set of 15 algebraic and linear partial differential 
equations to be solved. The number of equations may be reduced 
through the elimination of variables. There results a single 
diagnostic equation with geopotential height as the dependent 
variable. We develop a diagnostic equation for the geostrophic, 
hydrostatic adjustment first and then include the contribution from 
the radiative transfer equations. Two Lagrange multipliers are 
eliminated by combining (6), (7), and (8). Then, forming the 
vorticity from (3) and (4) and combining with (8) gives. 


7i 1 V 2 <J)-7i 3 <|)+A 20 -7i 1 (v°-u°) +7t 3 4>°+n 1 (F 5x +F 6y ) =0 


( 11 ) 
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Reducing the thermodynamic variables is done as follows. 
Divide (9) by y and operate by a. Eliminate brightness temperature 
between (1) and (10). There results two equations, 

l[^(r-ni4 !t -£Vj-o (12) 


where , 
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and, 
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( 14 ) 


Combining (12) with (14) and substituting (2) gives, 


T y 2 f=i { V 

n. 


( 15 ) 


» 3 - 1 


Eliminating the Lagrangian multiplier between (11) and (15) yields 
a diagnostic equation in the geopotential height. 
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where , 


F-n 1 (v°-u y °) +* 3 <t>°+ (-^r°) 0 +g e (F 5x +F 6y ) 


(17) 


Much effort was spent programming for (16) . The resulting 
solution was not considered satisfactory. Given the complex 
coefficient structures and the delicate convergence criteria, much 
additional effort was expended through six subsequent derivations 
to express the variational formalism in forms easier to understand 
and easier to solve. These efforts eventually led away from a 
direct inclusion of the radiative transfer equation in a 
geostrophic, hydrostatic atmosphere. 

In retrospect, it seems that the solution could have been more 
easily obtained if the 15 equation set was reduced to the following 
3 equation set: 

7i 1 V 2 <t>-n 3 <| > --X 20 +n 1 (v°-u y ) (F 5jf +F 6y ) (18) 


T—- ^(4>„+P) 


( 19 ) 
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These equations can be solved iteratively by first setting X 2 =0 
and solving (18) for $. Then (19) is solved for T and the 
temperature substituted into (20) to derive an updated X 2 . Then 
the updated values are entered into (18) and the cycle repeated 
until a satisfactory level of convergence is attained. 


3 . Results 

There are three important points to consider regarding (18) - 

( 20 ) . 

a) NO RETRIEVAL OF TEMPERATURE IS REQUIRED TO BLEND SATELLITE 
OBSERVED BRIGHTNESS TEMPERATURES INTO THE METEOROLOGICAL DATA 
MAINSTREAM. This means that none of the problems associated 
with temperature retrievals will be encountered. The integral 
term appears on the right hand side of (20) not as a term to 
be solved. This is analogous to solving the radiative 
transfer equation for the brightness temperature - a very easy 
exercise. This single finding may make it worth while to 
pursue the formal variational approach to assimilation of 
microwave channel data especially if higher resolution 
radiance data becomes available in the future. 
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b) There must be observations of geopotential height or winds or 
both of equivalent accuracy with the satellite measurements in 
order for the (18) - (20) to work. Accurate observations of 
temperature apart from space-based measurements are not 
necessary. The caveat is that geopotential height must be 
known at the boundaries of the domain in order to obtain a 
solution for (18) . Lateral boundaries would vanish for the 
equations written on the sphere and the top boundary 
conditions can be removed to the top of the atmosphere or to 
some level where model predictions or climatology give 
satisfactory estimates. 

c) It is highly probable that (18) - (20) converge to a solution. 
The same equation set with the absence of the second term of 
(20) (the radiative transfer equation contribution) is known 
to converge. The second term of (20) is an integral term 
which should further stabilize the solution. 

The main goal of the variational assimilation project was to 
blend satellite-derived thermodynamic data into the meteorological 
data mainstream in a dynamically consistent way. The classical 
variational calculus method used to achieve that goal typically 
yields sets of complicated equations that require innovative 
methods for solution and also involve immense programming efforts. 
Therefore the effort was broken down into several simpler models 
that could be solved. 
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The attempt to reduce the equation set from 15 equations to 
one diagnostic equation in geopotential height resulted in equation 
(16) . After an extensive programming effort, a satisfactory 
solution was not obtained. I was unable to devise a scheme that 
could determine whether the problems were mathematical or 
pr ©grammatical. During the six other efforts to derive a more 

tractable diagnostic equation a conceptual error was made, namely, 
was treated as a variable that could be differentiated with 
respect to a. The observed brightness temperature dropped out of 
the equations. This led to the conclusion that the satellite data 
could not be successfully included in a classical variational 
assimilation. 

With the discovery of this error during the writing of Chapter 
5 of this final report, that conclusion is no longer valid. It 
appears, instead, that the satellite data can be successfully 
incorporated into a variational assimilation and that the blending 
can be done without any of the problems typically encountered with 
temperature retrievals. 



